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ABSTRACT 

Arcsecond-resolution spectral observations toward the protobinary system 
IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillime- 
ter Array Complex organic molecules such as CH3OH and HCOOCH3 were 
detected. Together with the rich organic inventory revealed, it clearly indicates 
the existence of two, rather than one, compact hot molecular cores (< 400 AU in 
radius) associated with each of the protobinary components identified by their 
dust continuum emission in the inner star-forming core. 

Subject headings: astrochemistry — ISM: abundances — ISM: individual (IRAS 
16293-2422) — ISM: molecules — radio lines: ISM — stars: fomation 



1. Introduction 

The low-mass protostellar source IRAS 16293-2422 (hereafter 116293) is located in the 
p Ophiuchus cloud complex at a distance 160 pc from the Sun. High angular-resolution 
observations revealed that 116293 is also a protobinary system of two components with 
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a projected separation of ~5.2" (~840 AU), and 3-mm continuum observations (Looney, 
Mundy & Welch 2000) suggest that the northwest component (0.61 M Q , hereafter I16293B) 
is slightly more massive than the southeast component (0.49 M , hereafter I16293A). Single- 
dish submillimeter line-surveys led van Dishoeck et al. (1995) to conclude that, within a 
20" region, 116293 consists of a cold (Tk m — 10—20 K) outer molecular envelope, a warmer 
(Tkin — 40 K) circumbinary envelope of 10"— 15" in size, and a hot (Tkm > 80 K) region of 
dense gas only 3"— 10" in size. This innermost gas was found to be rich in organic molecules 
and it was suggested that their presence was due to outflows (Hirano et al. 2001) interacting 
with the circumbinary envelope. However, Ceccarelli et al. (1999, 2000) showed that the 
emitting region is warm (100 K), dense n H2 > 10 7 cm -3 , and very compact (~2"— 3"). This 
led to the suggestion that 116293 contained a hot molecular core about 150 AU in size where 
high dust temperatures (Td us t ^ 100 K) can facilitate the evaporation of icy grain mantles 
(Ceccarelli et al. 2000). This conclusion is supported by recent observations which show that 
116293 inner hot core has a molecular inventory rich in complex molecules, similar to the 
hot molecular cores (HMCs) associated with massive protostars (Cazaux et al. 2003). 

Nevertheless, several questions remain that cannot easily be answered by single-dish 
observations. One major concern is that its inferred size of ~150 AU (< 1") (Ceccarelli 
et al. 2000; Schoier et al. 2002) is much smaller than the binary separation in 116293, and 
comparable to the dimensions of protostellar disks. Furthermore, it is unknown whether 
high resolution observations of low-mass hot cores will exhibit clear chemical differentiation 
on small spatial scales, similar to that seen in massive HMCs (e.g. Wyrowski et al. 1997). 

Submillimeter Array 6 observations with arcsecond resolution (1"=160 AU) were carried 
out on 2003 March 14 (compact configuration) and July 12 (extended configuration) with 
5 antennae. The phase tracking center of the observations was a(J2000) = 16 /l 32 m 22 s .91, 
5(J2000) = -24°28'35".52. The digital correlator at the time was configured with eight over- 
lapping segments ("chunks") each of 104 MHz bandwidth and 128 channels, except one 
chunk which had 512 channels, giving a frequency resolution of 0.812 and 0.203 MHz, re- 
spectively. The total bandwidth of each sideband was ~670 MHz covering two frequency 
ranges: 343.555-344.225 GHz (LSB) and 354.211-354.881 GHz (USB). Quasars 1743-038 
and NRAO530 were observed for phase and amplitude calibration. The flux scale is esti- 
mated to be accurate to 25%. The synthesized beam sizes are ~l".3x2".7 at 344 GHz and 
~l".lx2".5 at 354 GHz with natural weighting; these beams correspond to a linear scale of 
~200x400 AU at a distance of 160 pc. The data were calibrated using the MIR software, 
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and imaging was done with MIRIAD software. 

2. Results and Discussions 

To make positive identifications of the spectral features detected, in addition to the 
basic requirement of frequency coincidence, an iterative process with the following stringent 
selection criteria was applied to each candidate species: 1). Expected line intensity - unfavor- 
ably weak transitions were rejected. 2). Looking for transitions with similar or stronger line 
intensities in the bands - if missing, the particular candidate line was excluded. 3). Looking 
for slightly less favored transitions of the same molecule for a possible sign of low-level (< 
3-<t) "detections" if no other stronger transitions predicted were available in the bands. This 
would be an indirect verification of the particular line identified. 4). Literature check against 
lists of known molecules in 116293 in the literature - not mandatory for a detection as SMA 
is sensitive to dense hot cores. 5). Making the requirement that the molecular candidate 
has to be a known hot core molecule. 

Table 1 lists some of the representative molecular transitions detected. The LSR velocity 
of each molecular line measured is also listed in Table 1. Single-dish observations reveal that 
the LSR velocity of 116293 is centered at 3.9 to 4.2 km s -1 and is often found within the 
range of 3.3—4.7 km s -1 . Asymmetric double-peak line profiles with the blueshifted peak at 
2.5—3.5 km s _1 and the red one at 5.0—6.0 km s _1 are readily seen in single-dish observations. 

On the other hand, an interferometer like the SMA can easily filter out the diffuse 
ambient cloud, such as the warm gas envelope of 116293 at V L sr ~3.9— 4.2 km s" 1 , and its 
small beam will tend to pick up the fine kinematic details of dense hot regions. The SMA 
beam is small enough to resolve the two protobinary components in 116293, thus SMA is 
able to trace individual velocity component in the rotating disks and/or hot cores. As shown 
in Fig. la, the two protobinary components are clearly resolved in the continuum emission. 
Sample images of five organic molecules are shown in Figs. lb-f. Fig. 2 shows the sample 
spectra of four large organic molecules plus HC 15 N in I16293A and I16293B. By examining 
HC 15 N emission toward I16293A, two major velocity components, one at 1.7 km s _1 and the 
other at 4.7 km s" 1 were seen (Fig. 2a); in I16293B, one component at 2.5 and the other at 
6.0 km s _1 were also exhibited (Fig. 26). Furthermore, the small beam of an array suffers 
less beam dilution, thus is sensitive to high-velocity components within compact sources. 
In the case of HCN 4—3 and HC 15 N, for example, linewidths >8 km s" 1 are seen toward 
I16293A. The HCOOCH 3 line in I16293A is also fairly wide, with a linewidth of ~8 km s" 1 
(Fig. 2e). The large linewidth (5—8 km s _1 ) is common for simple molecules like HCN, or 
strong molecular emission like HC 3 N and CH 3 OH. For weaker emission, typical linewidths 
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are of 2—4 km s 1 in I16293A; however, toward I16293B, the linewidths of most molecular 
lines are notably narrower, i.e., between 2 and 4 km s _1 . 

Because of array characteristics and sensitivity, it is not unusual that an array can 
preferentially pick up only the stronger emission peak of a double-peak profile as seen by a 
single-dish. The CH 3 OH line shows Vlsr ~ 2.5 km s" 1 toward I16293A, with two velocity 
components visible near 4.6 and 1.8 km s _1 (Fig. 2c). The kinematic structure traced is 
thus similar to what seen by HC 15 N; a similar line profile of 13 CH 3 OH was also observed 
toward I16293A. Toward I16293B, the CH3OH emission is weaker and is mainly from the 
blueshifted component near 2.5 km s _1 (Fig. 2d), which again agrees well with what was 
seen by HC 15 N. Since the optically thin HC 15 N emission basically traced the densest regions 
such as the disks, it suggests methanol may also reside in the same dense regions. However, 
it is also plausible that the two velocity components seen in the methanol line profile are 
due to self absorption at ~3.9 km s -1 . Toward I16293B, the LSR velocities of HCOOCH3 
lines were found consistent with each other at ~3.9 km s _1 , within the velocity uncertainty 
of ~0.7 km s _1 . However, a non-trivial velocity scattering with Vlsr between 3.9 and 6.0 
km s _1 was observed toward I16293A. A close examination of the integrated intensity maps 
of HCOOCH3 show their peak emission positions are situated at ~0".5 east (Vlsr — 3.9 km 
s _1 ) and south (Vlsr — 5.9 and 6.0 km s -1 ) of the I16293A center position where Vlsr ~ 
4.6 km s -1 . Again, such a velocity field variation traces the general velocity gradient seen 
by the HC 15 N emission nicely. 

All spectral emission appears to originate from two compact regions of a beam-convolved 
size of ~200x400 AU in radius, concentrated toward the two protobinary components 
I16293A and I16293B. An exception is the 4—3 HCN transition which shows complicated 
kinematic structure (Takakuwa et al. 2004, in preparation). 

The beam-averaged line-of-sight column densities and fractional abundances of all ob- 
served molecules are listed in Table 2, where a gas column density Nu 2 = 1.6 x 10 24 cm -2 
(Schoier et al. 2002) was adopted for both protostellar cores. Column densities were measured 
at the peak emission positions of each integrated spectral line and were derived assuming 
thermalized level populations and optically thin lines. At these high densities, gas and dust 
are expected to be thermally well-coupled and so an excitation temperature T cx ~ T dust = 
100 K was adopted, except for molecular transitions with energy levels > 500 cm -1 (>700 
K) where T ex = 300 K was used (Schoier et al. 2002). For HCOOCH3, with multi-line de- 
tections, the column densities were derived from rotation diagrams (Fig. 3); we found T cx 
= 105 K for I16293A, and 116 K for I16293B. The assumption of T cx = 100 K applied for 
transitions of lower energy levels is hence warranted. 

The derived fractional abundances and relative column densities are in good agreement 
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with those derived from submillimeter observations of the Orion KL and Sgr B2 massive 
HMCs (Sutton et al. 1991, 1995). When molecular emission is mainly from the compact 
cores, abundances derived using the smaller SMA beam can be much higher (e.g. 0C3H2) 
than found in previous single-dish observations (van Dishoeck et al. 1995; Cazaux et al. 2003). 
This explains why CH3OH abundances derived from modeling higher energy transitions 
(Schoier et al. 2002) in the inner 116293 core show good agreement with our values, and 
suggests that the lower array abundances for HCN and HCOOCH 3 is indicative of less 
compact emission where the extended component is filtered out by the interferometer. 

The high abundances of organic molecules, particularly methanol, indicates that icy 
mantles have recently been evaporated (Charnley, Tielens & Millar 1992). For many 
molecules detected in both sources, the measured abundances only differ by factors of about 
2—3; as expected if both cores collapsed from the same cold core material. Nevertheless, 
the apparent absence of c-C 3 H 2 and CH 2 CO in I16293A, if confirmed, raises the prospect 
that future observations could discover pronounced chemical differentiation between the two 
cores. 

We have resolved the hot cores surrounding both protostellar sources in 116293 and 
measured the chemical composition in each. These preliminary observations demonstrate 
that the SMA will be an important tool to further explore the connection between the 
volatile organic chemical composition of the hot cores associated with massive protostars, 
those of Solar mass, and the composition of comets. 

This work was supported by NSC 92-2112-M-003-006 grant (Y.-J.K.) and by NASA's 
Exobiology and Origins of Solar Systems Programs through NASA Ames Cooperative Agree- 
ment No. NCC2-1412 (S.B.C.). 
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Fig. 1. — Spectral images of large organic molecules toward 116293. (a) Continuum at 354 
GHz; (b) spectral emission of CH 2 CO; (c) c-C 3 H 2 emission; (d) CH 3 OH; (e) HCOOCH 3 at 
344029 MHz; and (f) CH 2 CHCN. Crosses mark the positions of I16293A and I16293B hot 
cores. The angular size for a linear scale of 400 AU is shown in (a). The grey scale denotes 
the 354 GHz continuum. Continuum emission was imaged with the wide continuum channel 
at 354 GHz, at a noise level of 0.11 Jy beam -1 , and the estimated line contamination is < 
3% in I16293A and < 1% in I16293B. The continuum flux density of I16293A is 4.97±0.49 
Jy, which is slightly lower than that of I16293B (5.14±0.61 Jy). The dark ellipse represents 
the HPBW of the synthesized beam. Contours are shown at 3-, 4-, 5-, 7-cr levels in general, 
then at irregular intervals up to the peak values; dashed lines indicate contours at negative 
3- and A-a levels. 
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Fig. 2. — Sample spectra of large organic molecules toward 116293. The left column shows 
spectra taken at I16293A; the right column, at I16293B. HC 15 N spectra are shown in (a) 
and (b); CH 3 OH, in (c) and (d); and HCOOCH3 (344029 MHz), in (e) and (f). (g) gives 
C 2 H 3 CN line and (h), CH 2 CO. Tentatively detected (CH 3 ) 2 1 72.i6-16i.i5 EA, AE, EE & 
AA transitions at 343753.3, 343754.2 and 343755.1 MHz can also be seen in (g). The weaker 
spectral appearance of the tentative (CH 3 ) 2 line is because spectrum (g) is taken at the 
peak position of the integrated CH 2 CHCN line. All spectra were Hanning smoothed for 
better S/N ratios. 
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Fig. 3.— The rotation diagrams of HCOOCH 3 . Note that the 344170 MHz doublet (Eu = 
454 K) appears to be situated at a location slightly higher in the rotation-diagram fits. It 
might imply this doublet is just in the transitional regime toward a higher excitation, or it 
might indicate simply this is a misidentified line. Excluding this doublet would result in a 
lower excitation temperature, which agrees better with the value (~60 K) of Cazaux et al. 
(2003). 
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Table 1. Sample molecules detected toward IRAS 16293-2422 protostellar cores. 



Molecule 


Transition 


Frequency 
(MHz) 


'-'Low 

(cm -1 ) 


Ii/,I16293A a 
(Jy bm- 1 ) 


V LS R,A b 
(km s~ 1 ) 


Ii/,I16293B a 
(Jy bm~ 1 ) 


V LS R,B b 
(km s~ 1 ) 


HCN 


4-3 


354505.5 


17.7 


10.88±0.44 


4.6 


9.92±0.44 


3.6 


HC 15 N 


4-3 


344200.3 


17.2 


5.26±0.24 


3.9 


1.72±0.24 


3.9 


C-C3H2 


23i3,io-23i2,ii 


343804.9 


548.9 


2-a c 




3.48±0.32 


4.6 


CH2CO 


172,15-162,14 


343693.9 


127.9 


2-CT 




1.98±0.28 


3.5 


HC3N 


39-38 


354697.4 


224.8 


3.30±0.20 


4.6 


1.54±0.20 


3.9 


CH3OH 


182,16-173,14 E 


344109.1 


280.0 


2.20±0.15 


2.5 


1.94±0.15 


2.5 


13 CH 3 OH 


4l,3-3 ,3 A 


354445.9 


18.6 


2.27±0.20 d 


2.5 


2-CT 




C2H3CN (Vl 6 =l) 


364,32-354,31 


343761.8 


564.9 


2.82±0.26 


5.3 


0.94±0.26 


5.3 


HCOOCH3 


32*, 32-31*, 31 A 


344029.6 6 


180.4 


2.28±0.16 


4.6 


2.96±0.16 


3.2 




322,3i-31i,3i E 


344029.8 f 


180.4 












28i8,io-27i8,9 A 


344170.9 


304.2 


1.09±0.20 


6.0 


1.04±0.20 


3.9 




28i8,n-27i8,io A 


344170.9 


304.2 












33*,33-32*,32 A 


354608.0° 


191.9 


1.71±0.20 


5.9 


2.75±0.20 


3.9 




332,32-32l,32 E 


354608.4s 


191.9 












128,5-117,5 E 


354742.4 


49.7 


1.85±0.19 


3.9 


0.98±0.19 


3.9 



a The line intensity measured in Hanning-smoothcd spectrum at the peak-emission position of integrated intensity map. 
This value is generally lower than the actual peak intensity shown in channel maps due to Hanning smooth. 

b The uncertainty of the LSR velocity determined is ~0.7 km s _1 . 

c Only a ~2-ct (< 0.6 Jy beam" 1 ) detection in channel maps. 

d Partially blended with weak HCOOH 6-type transition 17 ,i7-16i,i 6 at 354448.3 MHz. 
e * = 0, 1. 

f Together with transitions 32i,32-31i,3i E, 322,31-312, 30 E and 32i,32-3l2,30 E. 
g Together with transitions 33i,33-32i,32 E, 332,32-322,31 E and 33i, 33-322, 31 E. 
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Table 2. Molecular column densities and fractional abundances toward IRAS 16293-2422. 







I16293A 






I16293B 




HMC 


IRAS 16293 


Molecule 




N b 


X 




N b 


X 


X 


X 






(cm- 2 ) 


WJV H2 ) 




(cm- 2 ) 




WN H2 ) 


W^Vh 2 ) 


HCN 


91.80 


3.1(+14) 


2.0(-10) 


49.03 


1.7(+14) 


1.0(-10) 


3.2(-9) c 


1.9(-9) d 


HC 15 N 


42.86 


1.2(+14) 


7.4(-ll) 


6.97 


1.9(+13) 


1.2(-11) 




7.0(-12) d 


C-C3H2 








10.03 


7.2(+15) c 


4.5(-9) 


6.3(-ll) c 


3.5(-ll) d 


CH 2 CO 








3.22 


1.9(+15) 


1.2(-9) 


3.(-10) f 


1.8(-10) d , 5.0(-8)s 


HC3N 


16.73 


6.7(+14) 


4.2(-10) 


3.34 


1.3(+14) 


8.4(-ll) 


1.8(-9) f 


2.5(-ll) d , 1.0(-9)s 


CH3OH 


13.62 


1.1(+18) 


6.8(-7) 


6.20 


5.0(+17) 


3.1 (-7) 


1.4(-7) f 


4.4(-9) d , 3.0(-7) g 


13 CH 3 OH 


17.03 


8.1(+16) 


5.0(-8) 












CH2CHCN 


8.58 


1.5(+16) c 


9.4(-9) 


2.35 


4.1(+15) c 


2.6(-9) 


1.5(-9) f 




HCOOCH3 


h 


6.8(+15) 


4.3(-9) 


h 


4.2(+15) 


2.6(-9) 


1.4(-8) f 


2(-7Y 



a In Jy bm 1 km s 1 . 

b T cx = 100 K was assumed for all lines except where noted. For HCOOCH3, column densities were derived from actual 
rotation-diagram fits. 

c Sgr B2(M) hot core; from Sutton et al. (1991). 

d Single-dish observations with a ~ 20" beam; from van Dishoeck et al. (1995). 
e For transitions with energy levels > 500 cm- 1 , T cx = 300 K is adopted. 
1 Orion KL hot core; from Sutton et al. (1995). 
g IRAS 16293 hot core; from models of Schoier et al. (2002). 

h In I16293B, the integrated intensities of HCOOCH3 lines are 8.71, 3.14, 5.55 and 2.31 Jy bm -1 km s _1 with increasing 
frequency, and in I16293A, 16.22, 4.25, 5.12 and 4.22 Jy bm" 1 km s _1 . 

1 Single-dish observations with ~ 10" — 30" beam; from Cazaux et al. (2003). 



